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The theory of zonal polynomials is used to compute the average of a Schur polynomial of 
, argument AX, where A is a fixed matrix and X is from the real Ginibre ensemble. This 

C> generalizes a recent result of Sommers and Khorozhenko [J. Phys. A 42 (2009), 222002], 

and furthermore allows analogous results to be obtained for the complex and real quaternion 
Ginibre ensembles. As applications, the positive integer moments of the general variance 
Ginibre ensembles are computed in terms of generalized hypergeometric functions, these are 
written in terms of averages over matrices of the same size as the moment to give duality 
formulas, and the averages of the power sums of the eigenvalues are expressed as finite sums 
' of zonal polynomials. 

+2 . 1 Introduction 

a 

The statistical properties of the eigenvalues of random N x N matrices with independent, identically dis- 
tributed (i.i.d.) real entries is a prominent topic in both physics and mathematics. Physical applications 
began with the paper of May [231 on the stability of an ecological network consisting of many components 
' {2/j(^)}i=i....,JVi coupled in some unknown way. Suppose the evolution of the components is governed by 

a first order differential system, and write the linearization of the latter about a fixed point in the form 



(N 



O 



dyi{t) 



dt 



= (-lN+B)[yi(t)}i=i,...,N- (1-1) 

=1,...,N 



o 

With the N x N matrix B having all entries zero this differential equation exhibits exponential relaxation 
to the fixed point. For general B the system (|1.1[) is stable if and only if all the eigenvalues of B have 
real part less than or equal to 1. 

Since the coupling between components is unknown, it is reasonable to take the components of B to 
be random. May argued that for B a dilute matrix (fraction 1 — c of its elements zero) with non-zero 
elements i.i.d. random variables having mean zero and variance a 2 , the spectral radius will be less than 
1 provided o\J Nc < 1. This is an asympototic result, requiring that N be large. It is consistent with a 
rigorous result proved subsequently for the special case c = 1 [15] . 

In the mathematics literature, attention has focussed not only on the spectral radius, but also the 
eigenvalue density. With the elements standard Gaussians, by explicit calculation the eigenvalue density 
was proved to be asymptotically uniform in the disk of radius y/~N, centred about the origin in the 
complex plane [6] . This is the so called circle law [TTl |3] , which has recently been proved to remain true 
for general i.i.d. distributions [31]. Another mathematical property proved true in the Gaussian case 
by explicit calculation |7J is that for large N the expected number of real eigenvalues is asymptotically 
y/2N/it. Numerical evidence presented in [7] suggests this result persists for general i.i.d. distributions 
of mean zero and unit variance, although a proof is yet to be found. 
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Of interest in both the physics and mathematics literature has been the integrability properties of 
eigenvalue distribution in the Gaussian case. Ginibre |16j was the first to seek an analytic formula for 
the eigenvalue probability density function (p.d.f.), giving rise to the name real Ginibre ensemble for 
real Gaussian matrices. In [TB] an analytic expression was found for the eigenvalue p.d.f. conditioned so 
that all eigenvalues are real. For the conditioning specifying a general number of real eigenvalues, the 
analytic form of the eigenvalue p.d.f. was not obtained until the passing of a further twenty-five years 
[2"Tl [5] . And it has not been until the last few years that analytic computations based on the eigenvalue 
p.d.f. have been mastered to the extent that the probability of a prescribed number of real eigenvalues 
can be calculated [30] [TJ [55] , and closed form expressions for the correlations obtained [IT} [27J 03 HH1 HH] ■ 
Furthermore, these analytic studies have been extended (T2] [2] to the case of partially symmetric real 
Gaussian matrices [21) . 

A very recent result [28] relates to the average of the Schur polynomials s K (Xi, ■ ■ ■ , \n) with respect to 
the eigenvalues {Xj} of matrices from the real Ginibre ensemble. We recall that the Schur polynomials are 
the basis for symmetric functions of {A.,}, labelled by a partition k\ > n-i > • • • > k,n > of non- negative 
integers, given by the ratio of determinants 

det[A" fc+JV - fc ] ?fc= i N 
s K (X u . . . , Ajy) = ! JM (1.2) 

aet[Aj \j,k=i,...,N 

With {Xj} being the eigenvalues of the matrix X, the Schur polynomials are often written s K (X). Making 
use of knowledge of the explicit form of the eigenvalue p.d.f., it is proved in [55] that for X a member of 
the real Ginibre ensemble 

f » T((N~n + Kn + l)/2) 

(«,(*)>* = j r((iV-n + l)/2 ' allK " 6Ven (1.3) 

[ 0, otherwise 

where |«| := J2j=i K j- 

It is the purpose of this paper to give a different viewpoint on this result. Explicitly, (j 1 . 3 j) will be 
deduced as a consequence of the theory of zonal polynomials. This viewpoint will allow for analogues of 
(|1.3[) to be given for the complex and real quaternion Ginibre ensemble. It further leads us to the evalu- 
ation of the moments of the characteristic polynomial in terms of generalized hypergeometric functions, 
in the case that the Gaussian matrices have a (matrix) distribution with a general variance. This in turn 
allows us to express the moments as different matrix integrals, in which the size of the matrix is equal to 
that of the moments, giving duality formulas. In the final section, again for the three Ginibre ensembles 
with a general variance, the averages of the power sums of the eigenvalues are expressed as finite sums 
of zonal polynomials. 



2 Zonal polynomials 



An alternative characterization of the Schur polynomials (|1.3[) is as the eigenfunctions of the differential 
operator 

f/ A 9 y i N ~ l ^x 9 i 2 v X]Xk ( 9 ±.\ (2D 

M^dA.J + a 2^^ 9Xj + a ^ \j-\k\d\* dX k ) [2A) 

j=l J j=l J l<j<k<N J J 
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in the case a = 1, with the structure 

s K (Ai, ■ ■ . , Ajy) = m K + ' s ^ j (2.2) 

/J<K 

In (|2.2[) m K denotes the monomial symmetric function indexed by the partition n (e.g. with N = 2, = 
A1A2 + A1A3 + A2A3), /i < k refers to the dominance ordering on partitions specifed by the requirement 
that y^j— 1 _■ K j G = 1) • • • t-^)) an( i the a Kfi are scalars. For general a the eigenfunctions of 

(|2.ip with the structure (|2.2p are the symmetric Jack polynomials P^^ a \\i, . . . , Ajv)- 
Associated with a partition k is the generalized Pochhammer symbol 

N 



(a) 



n r( ;r V^A"' - (2 - 3) 

7=1 r^-u- 1 )/") 



and in terms of this and the classical Pochhammer symbol (it) n :— u[u + 1) • • • (it + n — 1) one defines 
_ a^[(N-l)/a + li a) - 1/a _ ^ (l + (j-z-l)/a + ^-^) 1/a 

K_ F^) ' J ( ) '" 1 <i|< W (l + (j--i-D/«)l/a ■ ^ 

The quantity in turn is used to define the renormalized Jack polynomials 

•\ k \\k\\ 



rW(\ x ) _______ 

k (Ai, . . . , An) - — -j. — 



P^ ) (A 1) ...,A W ). (2.5) 



In the cases a = 2, 1 and 1/2 the renormalized Jack polynomials are the so called zonal polynomials 
associated with the symmetric spaces gl(N,R)/0(N), gl(N,C)/U(N) and u* (2N)/Sp(2N) [22] . 

For present purposes, a key propery of the zonal polynomials is their appearance on the right hand 
sides of the matrix integrals [T^l [HI [__j . 



d 2 \AA T ) 



A = 2k 



( S x(AO))oeo { N) = \ d 2) ((in' ( 2 - 6 ) 

0, otherwise 

(s x (AU)s K (UUi)) UeU{N) = 6 X , K " (2.7) 

C K ((1) ) 

•> A ij.vi> s . ,,,,,, = ; ci 1/2 \(i) N y a = k2 ' (2-8) 

0, otherwise 

where in (|2.6p the partition 2k is the partition obtained by doubling each part of k, while in (|2.8p . k 2 is 
the partition obtained by repeating each part of k twice. On the left hand sides the averages are over 
the classical groups O(N), U(N), Sp(2N) of unitary matrices with real, complex and real quaternion 
elements respectively, endowered with the corresponding Haar measure. 

We immediately observe a structual similarity between (|1.3|) and (]2.6p . In fact, as will be shown in 
the next section, (|2.6p implies and furthermore generalizes (jl.3[) . 
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3 Averages over the Ginibre ensembles 



That zonal polynomials are intimately related to averages over matrices with Gaussian entries is the theme 
of the monograph by Takemura [3D] . This theme has further been developed in the works [TH [25] . More 
generally, zonal polynomials can be related to any measure d[i(X) on the space of random matrices with 
real, complex or real quaternion entries possessing the property of being invariant under the mappings 
X i ^ UX, X i-> XU for U e O(N), U(N), Sp(2N) respectively. With {-) x denoting an average in such a 
setting, and (•)[/ denoting the average over the corresponding classical group, to apply this theory to the 
Ginibre ensembles we first make note of a fundamental factorization property of the zonal polynomials 
with respect to the former, namely (see e.g. [22] ) 

C K ((1) ) 

Here d a \(l) N ) := C^\x\, . . . , Xn)\ X i=— =x n =i an d {a%}, {h} are the eigenvalues of A, B respectively. 
Proposition 1. \30f One has 

{CM(AXBrf)) x = ^^^M {C^{XX^) X . (3.2) 
(G K ((1) )) 

Proof. For any /(X) integrable with respect to dfi(X), the invariance of dfi(X) under X i— » UX tells 
us that 

(f(AXBX^)) x = ((f(AUXBX^U^))u)x, (3.3) 
while the invariance under X i— > XU gives 

(f(AXBX^)) x = ((f(AXUBU^X^))u) x . (3.4) 

Choosing / = C [ K a) in (JO]) and using (13.11) shows 

(C<?\AXBXi)) x = °S^ N M a HXBXl)) x . (3 .5) 
Choosing / = ct ] in and 

again using (|3.ip allows the right hand side of (|3.5p to be evaluated, and 
results. □ 



We can use (|3.2I) combined with (|2.6[) — (|2.8[) to obtain a generalization of (|1.3|) for each of the real, 
complex and real Ginibre ensembles. First we consider the case of a general measure dfi(X) invariant 
under multiplication by unitary matrices. 

Proposition 2. Let dfj,(X) be as required for f3.2\) . For real matrices 

MAX))x = { §^ {C< ' m - (3.6) 
0, otherwise; 



for complex matrices 



( Sfi (AX)s K (X^A^)) x = S X , K C ^ AA ^ (cP(XXi)) x ; (3.7) 

(CP ((1) N )) 2 
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for real quaternion matrices 



{s^(AX)) x =\ (C£ /2) ((1)")) 2 ( 3 - 8 ) 
0, otherwise. 



Proof. Consider first (|3.6p . It follows from (|2.6p that 

MAX)) X = ((s^AXO))o)x 

{C^{A T AXX T )) X _ 2 

0, otherwise. 

Making use of (|3.2[) in the case a — 2 gives (|3.6p . The derivation of (|3.7p and (13 . 8|) begins with (|2.7|) and 
(|2.8p . then makes use of (|3.2p in an analogous manner. □ 

We see from p. 61) and (|2.3p that (| 1 . 3[) is reclaimed if we can show that for dfi(X) oc e-fifflTrrf^^ 

1 (C^{XX T )) X = 2\^[N/2]W. (3.9) 
d 2) ((l) w ) 

For this we change variables XX T = A, using the result (dX) oc (det A)~ 1 / 2 (dA) (see e.g. [HI eq. (3.30)]), 
to obtain 

(CW(XX T )) X = ± f e -(V2)^A (det A yift c (2) {A) {dA) (310) 

where C is such that the RHS equals unity when k = 0^, and A > denotes that the integral is over 
the space of positive definite matrices. Changing variables now to the eigenvalues and eigenvectors of A 
shows that (|3.10p is proportional to 

/ d\ x --- d\ N T\\; l/2 e- x ^C < ?\\ 1 ,...,\ N )T\\\ k -\ j \. (3.11) 
Jo J ° i=i j<k 

We recognize the integral in p. lip as appearing in the integration formula (see [9l eq. (12.152)]) 
- d\x-- d\ N l[\?e- x 'Cl a \\ ll ...,\ N )l[\\ k -\ 3 \ 2 / a 

= Ci a \(l) N )[a+(N-l)/a + l]^ (3.12) 

(take a = 2, a = —1/2 and change variables ti i— s- ti/2), which itself is a limiting case of a generalization of 
the Selberg integral conjectured by Macdonald and proved by Kadell and Kaneko (see [13] and references 
therein), and (|3.9p follows. 

Proceeding similarly, making use of (|3.12|) for a = 1, we can show that for X complex, with d/i(X) oc 
e-^&x'HdX), 

^^ ) ^ 1) ( XX '))x = [N^. (3.13) 

And for X real quaternion with d^{X) oc e - Tl '( xx1 ) (<jX) — the trace now being with respect to the 
quaternion structure and so selecting only one of the diagonal elements from each 2x2 block — we have 

M1/2L um (C { ' /2) (XX'))x = 2-l«l[2Al 2 ). (3.14) 
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Substituting l[3"30j) . ((37T3]) and (|3TT4|) in Proposition [2] gives the sought generalization of (| 1 . 3|) for the 
Ginibre ensembles. 

Corollary 1. With the distribution of the real Ginibre ensemble proportional to e"! 1 / 2 ) 11 ^ ^ an d ffe 
distribution of the complex and real quaternion ensembles proportional to e~ TT ( xx1 \ one has for the real 
Ginibre ensemble 



fe(«)).v = { C< 2 >((1)») " ' ' (3.15) 
0, otherwise; 

for the complex Ginibre ensemble 

(s,(AX)s K (X^A^) x = iff^ CP(AAi); (3.16) 



and /or the real quaternion Ginibre ensemble 

2-1*1 [2iV]! 1/2) 



(^(A*))* = ^ ci 1/2 \(l) N ) K V " " (3-17) 
0, otherwise. 



The fact that Prop. 12.23] 



p (a) ((1) iV ) = ^M# (3 . 18) 
h K 



where [SJ Prop. 12.28] l/h K is the coefficient of (xi + • • • + xjv)' k in P« (x) allows the results of Corollary 
Q] to be written 



( Sfl (AX)) x = 



h K P^(AA T ), fi = 2K 
0, otherwise; 

( s ^AX)s K (X^A^) x = S^ K h K P^ 2 \AA^) 

{s»{AX))x = < . 

I 0, otherwise. 

respectively. In this form Corollary Q] appears in the unpublished manuscript [25] of one of us (EMR). 

4 Hypergeometric functions 

We know from workings in [4] that 

2l K l[iV/2]L 2) _ 1 



1 



c«( ( m"N! ( kIq=iJ 

2-1*1 [27V] l 1/2) _ _J_ 

C (i/ 2 ) ((1)W) -2l«l| K |! rf * 2|Q = 1 ' ^ 
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which we substitute into the results of Corollary [T] as appropriate. To see the consequence of this, we 
recall the definition of the generalized hypergeometric functions 

p F& (a 1 ,...,a p ,b 1 ,...,b q ;x 1 ,..., x N ) := £ ± t\£)' ' '^M <% a) (x u ...,x N ). (4.2) 

k [&i]« 

For p = 0, q = 1 this is a generalization of the binomial expansion and we have [9j eq. (13.4)] 

N 

1 F^ ) (a- 1 x 1 ,...,x N ) = l[(l~x ] r a . (4.3) 
i=i 

It follows from this that upon multiplying both sides of (|3. 15[) and (|3.17[) by [—r]^ /d'\o,—i, and both 
sides of (|3.16[) by [— r]^ [— r]«^/d' | a =i^| =i) an d summing over jj, (or /i and k in the case of (|3.16[) ), 
the left hand sides of each of the identities can be summed according to (|4.3|) . For the resulting matrix 
averages to be well defined, we require r 6 Z>o- To simplify the right hand sides, we note from (12. 3|) that 

= 2 2 I«I[«/2]M[(« + l)/2]< 2 >, [«]g? = [«]^[«- 1]^, 

then make use of (|4.2[) . Consequently we obtain the following set of matrix integral evaluations. 

Corollary 2. Let r G Z>o. One /ias, /or the real, complex and real quaternion Ginibre ensembles 
respectively, 

(detQljv - AXy) x = 2 F (2) (-r/2, (-r + l)/2; 2AA T ) 
(dct(I N - AX) r det(I N - X^A^ r ) x = 2 F (1) (-r,-r;AA t ) 

(det(I 2N -AX) r ) x = 2 F (1/2) (-r,-r-l;^t/ 2 ). (4.4) 

Suppose A in (14. 4|) is invertible, and write E = (A^A) -1 . We can then rewrite (14. 4| to read 

(detQljv - xl) r )x = 2 ^o (2) (-r/2, (-r + l)/2; 2x 2 E) 
(det(Ijv - il) r det(Ijv - ^ f ) r )x = 2 ^ (1) (- r > - r 5 M 2s ) 

(det(I 2JV - xXy) x = 2 F (1/2) (-r, -r + 1; M 2 E/2). (4.5) 

Here the averages are over N x N matrices X with real, complex and real quaternion elements respectively 
having a distribution proportional to g-^C^"^ 2 1 )/ 2 j n the real case and e — Tr ( xx ' s 1 ) j n the complex 
and quaternion real cases. Introducing Y = X — (1/x)In and setting W — YY^ we then have that W is 
distributed as a non-central Wishart distribution (see e.g. [H]) in the first case, and its complex and real 
quaternion generalization in the other two cases. Furthermore, for r — 2s even in the real case, and for 
general non-negative integer r in the complex and real quaternion cases, the determinants in (|4.5[) can 
be written entirely in terms of W, and we obtain 

|x| 2JVs ((det W) s ) w = 2 F (2) (-s, -s + 1/2; 2a; 2 E) 
\x\ 2Nr ((det W) r )w = 2 F (1) (-r, -r; |x| 2 E) 

|xr((detM/r/ 2 ) w = 2 F Q (1/2) (-r,-r-l;|x| 2 E/2). (4.6) 
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The identities (|4.6p are noteworthy for the fact that in the real case a different generalized hyperge- 
omctric function evaluation is known |241 Th. 10.3.7], 



((detWY) w = (det^^^ ^^+^^ FPC-aiJV^i-S- 1 ^ 2 ) (4.7) 
where 

rn 

T m (u) := tt™^- 1 )/ 2 TJ T(u - (i - l)/2). 

This implies (after some minor simplification) the identity between generalized hypergeometric functions 

2 F (2) (- S ,- S + l/2;y) = [iV/2]$(detY)\if V^iVA-Y- 1 ). (4.8) 

Using the property of Jack polynomials [9l Exercises 12.1 q.l and 2] 

(detYyP<?\Y- 1 )=p£\Y) 

where k s := (s — kn, s — fCjv-i) ■ . ■ ,s — Ki), noting from (|2 .4|) that 

2N ^ _ [(jV + l)/2]j 2 ; 
2^1 d; [(iV + l)/2]l 2) ' 

and using the property of the generalized Pochhammcr symbol ()2.3|) 



[u]<g> = (-1)1- lii ^ ^ 

[(7V-l)/a- U +l- S ]L Q) 

(a consequence of the functional equation for the gamma function), (|4.8|) can be verified directly by 
comparing coefficients of P^ (Y) on both sides. 

5 Duality identities 

There are many matrix ensembles of N x TV matrices {X} for which (det(Ijv — xX) r )x can be expressed 
in terms of an average over dual matrix ensembles where the size of the matrices is r x r (see e.g. [9]). 
An example of relevance to the present study is an identity of Fyodorov and Khorozhenko [H] (see also 
[8]), which reads 

(\det(zI N -AU)\ 2 P)ueu(N) 

oc/ dt,... dt P H \LU ' I H \ tk - tj \\ (5.1) 
J0 Ja i=i v (; i<i<fe<P 

To write the right hand side as a matrix average, we require a result from random matrix theory [9j 
Exercises 3.6 q.3] giving that the matrix Y = X^B-^^B-^X, where X is an (N + p) x p standard 
complex Gaussian matrix, and A is a complex Wishart matrix B = b^b with b and (JV + p) x (JV + p) 
standard complex Gaussian matrix, has eigenvalue p.d.f. proportional to 

n (i +tl w+2p n i**-*ji a - 
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Noting too that 

p 

Yl det(\z\\ + UAA^) = det(|z| 2 Ijv P + Y ® AA r ) 
i=i 

we see (|5.ip can be rewritten as the matrix average duality 



(| det(zI N - AU)\ 2p ) UeU[N) = (det(|z| 2 I Wp + Y <g> AA*)) Y . (5.2) 

The close relationship seen in Sj3] between averages over the unitary group which contain an arbi- 
trary matrix, and averages over the complex Ginibre ensemble, suggests analogous formulas hold for 
the averages in (|4.5[) . This is indeed the case. The sought identities follow from a representation for 



2-FJ) ( — r ' ~ a / a — ( r — 1))^")) r *= ^ + as an ^-dimensional integral, valid for all a > 0. 
Proposition 3. Let r E Z + , a > and Y be an N x N matrix. We have 
2 F^\-r,-a/a-(r-l);Y) 

-| pOO pOO r 

-. / dt 1 ---\ dt r l[e- t 'tf- 1 det(l N + (ti/a)Y) J[ \t k -tA 2a , (5.3) 

o-l,2a,r ■'0 JO l=1 Kj<k<r 



where 



/'OC /> oc n 

/ / dtn TT e-*'tf 1 ]"[ \t k -t 

Jo Jo ,_, ,^.^,.^„ 





■J I 



Z=l l<j<fe<n 

(^/ms normalization is a well known limiting case of the Selberg integral, and as such can be evaluated as 
a product of gamma functions Prop. 1^.1.3], although we don't need this fact). 

Proof. With k' denoting the conjugate partition, obtained by interchanging the rows and columns of 
the diagram of k, we have [9l Exercises 12.4 q.2] 

M W = ( _ a) -|«l[_ H (i/«). 

Also, from the definition of d' K in terms of arm and leg lengths [HI Eq. (12.37)], and the corresponding 
definition of h K [HI Eq. (12.58)], one sees 

d' K , = a |K| /i K U^i/a. 

Recalling (|2~5|) and ([472]) it follows 



2^o '(-r,-a/a- (r - 1); x u . . . , x N ) = ^ P^, ; (X). (5.4) 

""Kloii-tl/a 

The significance of this series form, as distinct from the form implied by (|4.2p . is that the coefficient of 
P^?\x) permits the integral representation 

a -\«\[ra]£ /a) [a+(r-l)a]£ /a) 
I a>— kL/qi 

i pOO pOO r 

— / dh-.- dtrHtf-'e-fp^iT) n \t k ~ tj \ 2a , 



Z=l l<j<fc<r 
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which is just a rewrite of (|3. 12|) . after making use of (|3. 18[) . Substituting in (|5.4p . and recalling the dual 
Cauchy identity for Jack polynomials [9, Eq. (12.186)] 

f[(l+ Xk y l )=J2P^(X)P! i } /a \Y), 

k,l=l K 

follows. □ 



Corollary 3. Let the averages over X be as in fr4.5\ ). We have 



(det(xI N - X) 2s ) 



x 



2 />OC />oc s 

/ dtx--- dt s \[ e~ u det(x 2 l N + i,E) [] (t k - tj) 4 

)A.r JO JO ,_ n i^„^i,^„ 



Wo,4,r JO 

(\det{xI N - X)\ 2r ) x 



1=1 \<j<k<s 



{dct{xl 2N - Xf) x 



-. / dtx--- I dt r Y[e tl det(|x| 2 l7v + i;E) JJ {t k -t 

0,2, r -/0 JO ; =1 K ,<Kr 



5= / dt x --- I dt r T\e- tl Aet{\x\H 2N + tiT) j~[ |t fc - i 

Wo,l,r -/O ./o , =1 X<j<k<r 



(5.5) 



(5.6) 



(5.7) 



Each of the right hand sides in the above identities can be written as matrix averages involving 
Wishart matrices. For example, with Y = a^a, where a is an r x r matrix of standard complex Gaussian 
entries, (|5.6[) can be written 



(| det(xI N - X)\ 2r ) x = (det{\x\% Nr + Y® E)) r . (5.8) 

In the theory of the complex Ginibre ensemble with E = Ijv it is well known that the eigenvalue 
density pn){{x, y)) (with TV i— > N + 1 for convenience) satisfies 

e-W 2 

= ^r<l det(zl w - *)| 2 )x (5.9) 

(see e.g. [TH eq. (1.3)]). Here z = x + iy and {X} is the complex Ginibre ensemble with a — 1^. Even 
though the equality is restricted to E = Ijv, this motivates us to use the duality identity (|5.6[) to investige 
the right hand side of (|5.9[) when {X} is the complex Ginibre ensemble for a more general variance matrix. 
In particular, suppose E = diag (cr, (l)^ -1 ). Then (|5.6[) gives that the right hand side of (|5.9[) equals 

For < |z| 2 < TV and iV large the leading behaviour is 

a + (1_^M 

7T 7T.V 

while for |z| = \/]V — r and cr fixed, for large N we obtain 

-!-(l +erf(V2r)) 
zir 
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independent of a. Since the boundary of the eigenvalue support is \z\ = y/N, this suggests our choice 
of S did not effect the largest eigenvalues. On the other hand, if we were to set a i— > ctVN, as well as 
\z\ = ^/N—r, the asymptotic form of (|5.10| is dependent on er, indicating that then the largest eigenvalues 
have been altered. 



6 Power sum averages 

In the previous section we computed, in (|4.5j) . the integer moments of the characteristic polynomials for 
real, complex and real quaternion Gaussian matrices with a general variance matrix E. Here we average 
the power sums Pk(X) := Tr (X k ) over the same class of matrices. For real Gaussian matrices with 
S = In, this average was computed using (|3 . 1 5[) in the case A = Ijy [28]. If we make use of each of the 
identities of Corollary [1] for general A, together with the expansion [22] Ex. 1.4.10] 

fc-i 

p fc (X) = ^(-l)' S(fc _ ui) (X) 

1=0 

we can generalize the result of [28] . 

Corollary 4. Let the averages over X be as in |^.5p . We have 

0, k odd 



fc-i [A r](i) 

0, fc odd or ft/2 > AT 

<»(*)>* =^ 2 - fc/2 C (l/2) (s) fc < 27Veven . 

^((iD 

Let and p°^((x, y)) denote the density of the real and (upper half plane) complex eigenvalues for 
real Gaussian matrices with variance matrix S. Then we have 

/OO P 
x k p\ l) (x) dx+ ((x + iy) k + (x- iy) k )p c (1) ((x, y)) alxaly 

where := {(x, y) : x £ R, y £ R + }. The average in Corollary [4] in the complex case can be rewritten 
as an average involving both the one and two point correlations, which we refrain from writing down. In 
the real quaternion case, the average again can be written in terms of just the one point density. Here 
there are no real eigenvalues, and we have 



(p k (X))x= / {{x + iy) k + {x-iy) k )p c {1) ((x,y)) alxaly. 
Jr 2 + 
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